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The fermion four-point functions and condensates as the chiral symmetry or- 
\ der parameters are calculated analytically in U(l) gauge theory in the massless 

phase. It is shown that in leading order of the loop expansion of the effective 
r S^ • action, there is a critical coupling for the nonperturbative parity-invariant 

chirality-changing four-fermion functions; this critical coupling however is the 
same as that for dynamical mass generation. Therefore there is no breaking 
r-j ! of chiral symmetry through the condensation of four-fermion functions below 

the critical coupling for dynamical mass generation for the fermions. 
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I. INTRODUCTION 



Dynamical chiral symmetry breaking (DCSB) plays a very important role in the research 
of fundamental interactions. Up to now, most investigations are focussed on DCSB with 
mass generation for the fermions. In such a scenario, the bilinear fermion condensate is taken 
as the order parameter of chiral symmetry breaking and it is known that it can receive a 
nonzero value when the coupling of the system exceeds some critical value. For a review on 
dynamical chiral symmetry breaking, see e.g. [|TJ. 

But this is not the complete story for DCSB, since the order parameter may not be 
the bilinear fermion condensate: in principle, higher n-point functions which break chiral 
symmetry might become nonzero, without dynamical mass generation for the fermions . 
In a massless U(l) gauge theory for example, four-fermion Green's functions which break 
chiral symmetry might become nonzero at values of the coupling below the critical value 
for dynamical mass generation. Such a breaking of the chiral symmetry by the four-fermion 
condensates, without mass generation for the fermions, was discussed recently by Holdom 
and Triantaphyllou Q, who have calculated these four-point functions numerically and found 
a critical coupling constant which is almost the same as that responsible for dynamical mass 
generation. 

One could also think of theories where some other symmetry prevents the dynamical 
generation of a bilinear condensate and a dynamical fermion mass, but where a nonzero 
condensate of a higher n-point function is not restricted by this symmetry. In a chiral gauge 
theory, the gauge symmetry will be broken by the generation of a fermion mass, but a 
chirality-changing four-fermion condensate does not break the gauge symmetry, and might 
therefore be generated dynamically. In such a scenario, the chiral symmetry can be realized 
in more delicate way which may be useful for model building j§-|[J: it allows for a hierarchy 
of chiral symmetry breaking, first by the four-fermion condensate at relatively small values 
of the coupling, and eventually a second phase transition, at higher values of the coupling, 
leading to mass generation for the fermions, due to the four-fermion condensates. For these 
purposes one should study this phenomenon in a chiral gauge theory, as suggested in || , but 
for simplicity and to see whether such a symmetry breaking patterns is possible, we restrict 
ourselves to a U(l) gauge theory. 

In our paper we confirm analytically the numerical results obtained by Holdom and 
Triantaphyllou @, using a slightly different method to calculate the nonperturbative four- 
fermion functions. We find that there is a critical coupling for dynamical chiral symmetry 
breaking via the four-fermion condensates, below which there are no nontrivial solutions for 
the four-fermion functions which break chiral symmetry, but preserve parity, Lorentz sym- 
metry, and the flavor symmetry. This critical coupling is (in our truncation scheme) exactly 
the same as the critical coupling for fermion mass generation. This result is obtained by 
calculating the loop expansion of the effective action to leading order for the nonperturbative 
four-point functions. Our result is also in agreement with an intuitive argument comparing 
the condensation of two fermions to the condensation of four fermions: for the four-fermion 
condensate we need at least the same binding force as for the bilinear fermion condensate. 

Our paper is organized as follows. In the next section, we discuss some general features 
about dynamical symmetry breaking, and we show how in principle the chiral symmetry can 
be broken without mass generation for the fermions. In Sec. |TJ, we derive a linear integral 
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equation for the four-fermion function, using the loop expansion of the effective action. Next, 
we analyze this equation for the nonperturbative chirality-changing four-fermion functions, 
using some expansions in order to perform the angular integrations analytically; more details 
can be found in the appendix. We find that there is a critical coupling for these functions; 
this critical coupling turns out to be the same as the critical coupling for dynamical mass 
generation. Finally we give some concluding remarks in Sec. [V[ 



II. GENERAL CONSIDERATIONS ON DYNAMICAL SYMMETRY BREAKING 

A. Symmetry Breaking Patterns 

The kind of symmetry breaking we are dealing with in this paper is the breaking of chiral 
symmetry in a U(l) gauge theory with Nf massless fermion flavors. The global continuous 
chiral symmetry is Ul(N/) <8> Uu(Nf), and this symmetry can be broken in different ways, 
characterized by the corresponding chiral order parameters. Since the flavor symmetry is 
not broken in vector-like gauge theories , only flavor-symmetry invariant order parameters 
may receive nonzero contributions; also parity || and Lorentz symmetry is believed not to be 
broken spontaneously, so in this paper we limit us to parity-, flavor-, and Lorentz-invariant 
order parameters. 

The simplest and most wellknown way of breaking the chiral symmetry is characterized 
by the trace of the two-point function (0 1 T-0 (xc) -0 (y ) 1 0) . This condensate, (0 |-0(cc)-0(a;) |0) , is 
the only possible two-fermion condensate which keeps the flavor symmetry and is parity and 
Lorentz invariant. It gives rise to a dynamically generated fermion mass, and the symmetry 
breaking pattern is 

U L {N f ) ® U R (N f ) — > U L+R (N f ) . (2.1) 

An explicit mass for the fermions would also break the symmetry in this way. 

However, this is not the only way to break the chiral symmetry, an alternative way is 
characterized by condensates of four fermions. The possible four-point function which keep 
flavor and Lorentz symmetry and which are parity invariant are 

Go,s+p(x 1 ,x 2 ,x 3 ,x 4 ) = (0\T4i(x 1 )\ ip(x 2 )4>(x 3 )\ 'ijj(x 4 )\0) 

+ (0\T^(x 1 ) l5 X°^(x 2 )^(x 3 h 5 X°^(x 4 )\0) , (2.2) 



Nf-l 



Ss+p{xi,x 2 ,x 3 ,x 4 ) = 



a=l 



(0|T^(x 1 )A^(x 2 )V(^)AXx 4 )|0) 



+(0\T$(x 1 )y s \ a i/>(x 2 )${x 3 )y s \ a i/>(x i )\0) 



(2.3) 



£o,t(zi,z 2 , x 3 , x 4 ) = (0|T^(x 1 )a^AV(x 2 )^(a;3)a^AV(x 4 )|0) , (2.4) 

Nj-l 

g T ( Xl ,x 2 ,x 3 ,x 4 ) = £ (0\T^(x 1 )a^X a ^(x 2 )^(x 3 )a^X a i>(x A )\0) . (2.5) 



a=l 



In these expressions, X a is the flavor matrix with a = l,...,Nj — 1; A is proportional 
to the unit matrix. These four-fermion functions can be related to the axial charge, and 
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the condensate of these four-fermion functions (with all the coordinates Xk in the same 
space-time point) can be considered as the order parameters of DCSB. 

A nonzero condensate of (one of) these four-fermion functions gives the symmetry break- 
ing pattern 

U L (N f ) (g) C^(JV/) — U L+R (N f ) ® Z, . (2.6) 

Here Z 4 is a nontrivial chiral subgroup which plays a role in preventing dynamical mass 
generation for the fermions. It corresponds to the discrete chiral transformation 

V>(ar) -> e^^^(x) , (2.7) 

and therefore the fermions remain massless in such a breaking of the chiral symmetry. 

Note that in principle there are a lot more four-fermion functions which break the chiral 
symmetry. In the case of one fermion flavor, Nf = 1, there are in total 54 independent 
functions which all lead to the same chiral symmetry breaking pattern, Eq. Q2.6| ), but the 
ones above, which we will call S + P and T type four-fermion functions, are the only ones 
which preserve parity and Lorentz symmetry. The other 52 independent functions are not all 
parity-invariant or Lorentz scalars, which means they will not receive nonzero condensates, 
assuming that the Lorentz symmetry and parity remains unbroken. However, it is not 
excluded that (some of) the nonlocal four-fermion functions become nonzero, in which case 
they might couple to the S + P and T type functions. We leave this problem to future 
investigations and for simplicity we only consider the nonperturbative four-point functions 
which are flavor-, Lorentz- and parity-invariant. 

This scenario can easily be generalized to condensates of 2n fermions, breaking the chiral 
symmetry into 

U L (N f ) <g> U R (N f ) — U L+R (N f ) <g> Z 2n , (2.8) 
with the nontrivial chiral subgroup Z 2n , corresponding to the discrete chiral transformation 

1>(x) -> e^ 75 ^(x) , (2.9) 

unbroken. 



B. Goldstone Bosons 

It is wellknown that if a continuous symmetry of the Lagrangian is broken dynamically, 
there will appear some massless bosons, the socalled Goldstone or Nambu-Goldstone bosons 



0,10]. These Goldstone bosons couple both to the conserved current which is associated with 
the symmetry, and to some fields which characterize the breaking of the symmetry. Note 
that these fields 0, which can be regarded as the order parameters of the phase transition, 
can be either composite or elementary. 

In the case of dynamical breaking of chiral symmetry, the Goldstone bosons are re- 
lated to the chiral-symmetry breaking n-point functions through the axial (or chiral) Ward- 
Takahashi identity. The conserved chiral current is 



4 



(2.10) 



and we consider the following Green's function in momentum space 

Gftq) = J d 4 xe^(O|jt(O)0(a;)|O) , (2.11) 

where is a (composite) field associated with the chiral order parameter we are considering^]. 
The breaking of chiral symmetry means that this order parameter becomes nonzero, and 
due to the chiral Ward-Takahashi identity this leads to a pole at zero momentum in the 
above Green's function; in the limit g M — > this Green's function behaves as 

G%(q) ~ q ^f dxe^id v (0\f 5 (0)^M (2-12) 
^Z|! (0 |^|0). (2.13) 

To be more specific, we first consider conventional chiral symmetry breaking, with a 
dynamically generated mass for the fermions. In that case, the composite field <fi(x — y) is 
i/j(x)^(y), and the chiral Ward-Takahashi identity leads to 

(p - k)»G$(p — k) — S(p) l5 + l5 S(k) , (2.14) 

where S(p) is the full fermion propagator. From this equation, we see immediately that a 
dynamically generated mass for the fermions will lead to a massless pole in the pseudoscalar 
Green's function G$(p — k) if p — > k; this pole corresponds to the massless Goldstone boson. 

In the case of dynamical chiral symmetry breaking through the four-point functions, 
as we are considering here, the situation is somewhat more complicated since we have to 
consider a composite field <fi which consists of four fermion fields. Using the chiral Ward- 
Takahashi identity, one can show that the Green's functions G 5 ^(q,pi) have to satisfy the 
relation 

%{G^(q,Pi)) abcd = (j5)aa>Ga>bcd(Pl + 1,P2,P3,P4) + Gab'cd(Pl,P2 + <?, P3, Pa) {l<j)b'b 

+ (l5)cc'Qabc'd(Pl,P2,P3 + Q,Pi) + Gabcd* (pi , P2 , Pz , Pi + q)(lb)dd' , (2.15) 

where Gabcdipi, P2, P3, Pa) is the Fourier transform of a four-fermion function, with momentum 
conservation Y^Pi — 0- In the limit g M — > 0, this leads to massless poles, if either the S + P 
type or the T type four-fermion function 

Gs+p(PuP2,P3,Pa) = J d 4 x 1 ---d 4 x 4 e- i ^ p ^(0\Ti ! (xMx 2 )i ! (x 3 )tP(x i )\0) 

+ J d 4 a; 1 ---d 4 x 4 e-^^(0|T^(x 1 )75V'(x2)V'(x3)75V'(x4)|0), (2.16) 
Gt(Pi,P2,P3,Pa) = J rf 4 x 1 ---rf 4 x 4 e- l Sp^(0|T^(x 1 )a^(^)^(^3)^(^4)|0), (2.17) 



1 Here we only write down the space-time coordinate x which is related to Goldstone pole in 
momentum space; in general, <j) can also depend on other space-time coordinates, which do not 
relate to the Goldstone pole. 
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is nonzero. This can be seen easily by decomposing the Green's function Gg into 

G$(q,Pi) = ^(xs+p(<?,ft) (I ® 7s + 75 ® /) + *r(g,Pi) (tb<v) ® o-^) + . . . (2.18) 
so Eq. (|2.15|) leads to 

xs+p(q,Pi) = Qs+p(pi + q,P2,P3,pi) + Gs+p(pi,P2 + q>Pz,Pi) 

+Gs+p(Pi,P2,P3 + q,Pi) + Gs+p(PuP2,P3,P4 + q) , (2.19) 

xr(q,Pi) = Gt{p\ + q,P2,P3,pi) + Gt(pi,P2 + q,P3,Pi) 

+Gt(pi,P2,P3 + q^Pi) + Gt(pi,P2,P3,P4 + q) ■ (2.20) 

Thus massless Goldstone bosons appear as soon as any of these nonperturbative four-point 
functions becomes nonzero. 

If we restrict ourselves to local composite fields <fi we get a similar expression, but now 
involving the condensates of these four-point functions. Defining the condensate as 

j^-p Gi(pi,P2,P3,P4)S(pi+p 2 +P3+P4) , (2.21) 

where i denotes either S+P or T, the chiral Ward-Takahashi identity relates this condensate 
to the Goldstone bosons. In momentum space this finally leads to massless poles for q 2 — ► 

Xs+p(0) = 4: J d Pl t^y 2 P3 Gs+p(pi,P2,P3, -Pi -P2-P3) , (2.22) 

Xr(0) = 4y — -— Gt(Pi,P2,P3,-Pi ~ P2 ~ Ps) ■ (2.23) 

Thus the condensates of the four-point functions we are considering here are directly related 
to the existence of massless Goldstone bosons. 

Note that both a nonzero condensate and a nonzero four-point function would signal 
dynamical chiral symmetry breaking, and in principle it is conceivable that the condensate is 
zero, and the four-point function not; the other way around (no nonzero four-point function, 
but only a nonzero condensate) is obviously impossible. In this paper, instead of four- 
fermion condensate, we directly calculate four-point function and take it as a generalized 
order parameter. 

The generalization to the case of N fermion flavors is in principle straightforward; also 
the generalization to higher n-point functions is a trivial exercise. 

III. EFFECTIVE ACTION 

In order to calculate the chirality-changing four-fermion functions, we need some non- 
perturbative calculation scheme. For this purpose, we use a loop expansion of the effective 
action |TTHl3|. For simplicity we will restrict ourselves to just one fermion flavor, so the 
generating functional in Minkowski space is 

Z[I, I, J] = e^^^l = J VfDfDA^ ^I^c(A)+m~e4^+i^i+J^] ? (3 
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where C(A) is the Lagrangian for the photon field 

C(A) = -\F, V F^ - £ {d,A»f , 



(3.2) 



in a general covariant gauge. In order to calculate the effective action, we start by integrating 
out the photon field, which gives 



J VA^ I ^^i^+i-^i^+J^M = 

- J d A xd A y (-e^(xh^(x) + r(x)) D^x, y) {-e^{y)Y^{y) + r(y)) 



cxp 



, (3.3) 



where D^x^y) is the bare photon propagator 



D^(x,y) = — 



5\x - y) , 



(3.4) 



in which the derivative is acting on the coordinate x, and a is the gauge-fixing parameter. 

Next, we introduce a generalized Green's function containing four space-time points, 
with spinor indices a, b, c, and d 



D ti{x, y, y') = d^ u (x, yha d jbc s \ x ' - y)^(y' - x) , 



(3.5) 



such that 



-D^ix, y) U(x)^(x)) U(y)i v ^(y) 

\ / aa \ 

= J d A x'd A y'D^(x, x', y, y')r{x)r{x'W{y)^ d {y') 



(3.6) 



Note that in the case with Nf fermion flavors, we just have to add additional indices to label 
the different flavors. Thus the generating functional can be written as 

Z[I, I, J] — J VipVi/j exp i J d A x (ijkfol) + ty + ipl) 

+ ^ J d A xd A x'd A yd A y'D^(x, x', y, yO^W V)^(vWV) 
~Jd A xd A yD^x,y) (-2e^(x)^(x)r(y) + J»(x)r(y 

Inserting a Gaussian integral of a bilocal auxiliary field Xab(x,x') 



(3.7) 



const 



J "^Xexp J d A xd 4 x'd A yd A y'Dlj(x,x',y,y') x 
( X ac (x,x') - er(x)r(x')) (x M M) ~ e^WV)) 



(3.8) 



gives for the generating functional 



7 



Z[I,I, J] 



V X DifjT>ifj exp i 



d 4 x(^l)(x)i$ x ip(x) + I(x)ip(x) + ip(x)I(x) 



+ I d*xd 4 x'd%d%'D^(x } x' } y, y') ( er(x) X bd (y, y'W{x) - ^ X ac (x, x') X bd (y, y') 



+e / dVy^x^D^x, y)r(y)^(x) - - / dVyJ^D^x, y)J v {y) 



(3.9) 



We see that the fermion field appears quadratically in the generating functional, so that 
they can be easily be integrated out. Thus we arrive at 



Z[I,I,J] = I Vxe iScAxJj ' J] 
T>x exp i ' 



dVx'd^y'D^x, x', y, y') X ac (x, x') X bd (y, y') 



+zTrlnS - / d 4 xd*y -J"(x)D. l/ (x,y)r(y) + I a (x)S ab (x,y)I b (y) 



, (3.10) 



in which S ab (x,y) is defined as 



^ac( X , x ) 



X + er J d^yD^x, y)J V {y))^ ^{x -x') + ej d 4 yd 4 y'D? d (x, x\ y, y') X bd (y, y') 
(i$ x + e^J d*yD, u (x,y)r(y)^j\x - x') + eD^(x, x'K d ^ bcX bd (x', x) . (3.11) 



The effective action for the %-field, induced by Eq. ( |3.10|) and with external sources, is then 
given by 

Se«[x, I, I, J] = ~yJ d A xd*x'D^(x, x'ha d lb cX ac (x, x') X bd (x', x) 

+zTrln5 - J d 4 xd 4 y (^(x)D fiV (x,y)r(y) + I a {x)S ab {x,y)l\y)^ , (3.12) 



where we have used Eq. (|3.5|). Given the effective action S e g[ X} I,I,J] in Eq. ( |3.12| ), the 
calculation of various physical quantities concerns the calculation of the following generating 
functional 



Z[I, I, J, J] = e 



iW[I,I,J,J] 



V X exp i 



S eS [ X ,I,I,J}+ d 4 X J B (X) XB (X) 



(3.13) 



where J B (X) is an external source coupled to the X b(X) field, and we use a shorthand 
notation X and B for the bilocal variables (x, x') and indices (ab). The generating functional 
Z[I, I, J] in Eq. Q is related to Z[I, I, J, J] by Z[I, I, J] = Z[T, I, J, J = 0]. Now we define 
the classical field X b{X) by 



Xb{X) = 



5W[I,I,J,J\ 
5J B (X) 



(3.14) 



S 



and the partial Legendre transformation 

T[I,I,J,x] = W[I,I,J,J] 
It is easy to see that 

sf[I,i,J,x} 



5xb{X) 



d^X J»(X)x B {X) 



J»{X) 



and hence 



W[I,I,J] = W[I,I,J,0] = f[I,I,J, X ] 



(3.15) 



(3.16) 



(3.17) 



Therefore what we need to do is to calculate T[I, I, J, x}- When we turn off the source for 
the bilocal field, J B (X), Eq. ( gTTBp becomes 







(3.18) 



which reflects the stability of the physical vacuum. The standard loop expansion formula 
for r[J, /, J, x) is given in |12|], and reduces in the present case to 



f[I,I,J, X } = S eS [x,I,I,J] + ~Trln 



s 2 s eS [xJJ,J] 

^Xsi (Xi)5xb 2 (X 2 



+ r 



2 - 



(3.19) 



where the first term is the tree level contribution, the second term is the one-loop contri- 
bution, and T 2 represent higher-loop contributions, containing the sum over all higher loop 
1PI vacuum diagrams with the propagator 



and vertices 



6 2 S cS [ X ,ll, J] 
8xb 1 (X 1 )8xb 2 (X 2 ) 

8 n S cS [xJ,I, J] 



8 X B x {X 1 )...5 X B n {Xn) 



(n > 3) . 



(3.20) 



(3.21) 



This expansion, Eq. ( |3.19|) , gives systematic rules to calculate the effective action, and 
we base our further discussion on this expansion. For simplicity, we consider the effective 
potential up to tree-level order only in our present paper, but in principle it is straightforward 
to go beyond this leading order calculation. At tree level, when the external source is 
switched off, we find 



W[I,I,J} = S cS [x,ll,J] 
at the stationary point of the effective action 

6S eS [xJJ, J] 



5x ac (x,y) 



0. 



(3.22) 



(3.23) 
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Using the definitions Eqs. ( |3.12 ) and (|3.11 ), the stationary condition for the effective action, 
Eq. ( ggp , leads to 

(3.24) 



X ac (x, x') = -ieS ca (x\ x)-ej d 4 yd 4 y' S cb (x' , y)I b (y)I d (y')S da (y' , x) . 
Combined with Eq. ( |3.11| ), we find 

S~ c \x,x') = ^ x ) ac 5\x-x') + e^ ac 5\x-x') I d A yD^x,y).r{y) 



tS db (x, x') + / d*yd 4 y'S da ' (x, y')I a ' (y')I c ' (y)S c ' b (y, x') 



j4„, j4„,/ ada! 



c'b, 



(3.25) 



Note from Eq. ( |3.22| ) and Eqs. ( |3.12| ) and ( |3.23j ), that the function 5" is related to the full 
fermion propagator (two-point function) 



8I a (x)5I b (y) 



Sab{x,y) 



d zlJz 



5S cb (z,y) 
5I a (x) 

Sab(x,y) + I d A z ^ ' I c (z) 



(3.26) 
(3.27) 



where the functional derivatives 6S/5I and 5S/5I can be calculated using Eq. ( p.25[ ). Com- 
bined with Eq. ( |3.24j ), we find that the x _ field is the full fermion propagator if we switch 
off the external sources. On the vacuum, all external sources should be taken to zero, and 
therefore Eq. ( |3.25| ) gives 



Sj{x,x') - (i$ x )a C 5 A (x - x') = -ie 2 D flv (x,x')^S dh (x,x')^ 



(3.28) 



which is just the quenched ladder approximation for the Schwinger-Dyson equation for the 
full fermion propagator. By differentiating the generating functional twice with respect to 
the source J and imposing the stationary condition for the effective action, we can get the 
Schwinger-Dyson for the photon field. 

For the four-point function, we use Eq. ( |3.26| ) and ( |3.27| ) to relate the four-point function 
to derivatives of the function S. After putting the sources equal to zero, and defining 



G acbd (x,x',y,y') 



SS ac (x, x') 



5I h {y)8I d {y>) 



j=o,/=/=o 



(3.29) 



we find for the full four-point function 



g acbd (x,x f ,y,y>) = 



8I a (x)5I c (x')5I b (y)5I d (y>) 



J=0,/=7=0 



G acbd (x,x',y,y')-G bcad (y,x>,x,y') 



bead I 



G 



acbd. 



x,x ,y,y 



G 



adbc. 



x, y ; y t x 



(3.30) 



We have written down all possible expressions, which we can use to solve the four-point 
functions. Note that the four-point function Q is written as the sum of the s- and t-channel 
components G. 

Using Eq. (|3.25| ) we find for the s-channel 
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d 4 yd 4 y'S; b \x, y)G bda ' c \y, y', x 1 ,x' 1 )Sg 1 (y', x') 



A* 

''lad 



iG 



dba'c' 



e 2 D 



fj, u {x, x J7 6c 



(3.31) 



This is nothing but the (inhomogeneous) Bethe-Salpeter equation. For the S + P and T type 
four-point functions in the massless-fermion phase (or more general: for all nonperturbative 
four-point functions), the second term on the right hand side of Eq. (|3.31|) vanishes, so we 
arrive at following equation, namely 



G a c (xi, X2, x 3} X4) — 

-ie 2 / d 4 yd 4 y'S ad/ (x 1 ,y)^, a ,G a ' b ' cd (y, y', x 3 , xM c ,S c > b (y', x 2 )D flv {y, y') . 



(3.32) 



This equation is the s-channel equation for the nonperturbative four-fermion functions in the 
quenched ladder approximation. For the i-channel equation, we can get a similar equation 
in exactly the same way. This t-channel equation however is completely equivalent to the 
s-channel equation, and it can also be obtained by applying a Fierz transformation to the 
s-channel equation. 

Next we make a Fourier transformation to momentum space. The equation for the 
s-channel four-point function becomes 

yabcd 



G abcd ( Pu p 2} p 3 ,p A ) 



ie 



I -^ S ^'(PM'a'G a ' b ' cd (pi + q,P2~ q^PiHM-p^D^q) , (3.33) 



(2*) 



where D fiv {q) is the bare photon propagator in momentum space 



D^{q) = 



'l + a 



%qv 

q 2 



(3.34) 



As mentioned before, the t-channel equation is completely equivalent to this s-channel 
equation. For the complete nonperturbative four-point functions in momentum space, 
Q(Pi,P2,P3,Pi), we thus find 



G abcd (PuP2,P3,P4) = G abcd ( Pl ,p 2 ,p 3 ,p 4 ) - G cbad (p 3 ,p 2 , Pl ,p 4 ) , 
= G abcd ( Pl ,p 2 ,p 3 ,p 4 ) - G adcb ( Pl ,p 4 ,p 3 ,p 2 ) . 



(3.35) 
(3.36) 



So finally we have the leading order equations for the nonperturbative four-fermion functions, 
which is shown in diagrams Fig. 0. For the perturbative four-fermion functions ones we have 
to add the tree contribution coming from the second term in the right hand side of Eq. ( |3.31| ) . 



IV. CALCULATION OF THE FOUR-POINT FUNCTIONS 

One can show that the perturbative four-fermion functions do not contribute at all to 
the nonperturbative ones, so the nonperturbative four-point functions form a closed set of 
equations, at least below the critical coupling for dynamical mass generation. Furthermore, 
since we assume that all the other types of nonperturbative four-point functions (e.g. the 
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ones which break parity or Lorentz invariance) remain zero, and we only have to deal with 
the S + P and T type four-point functions 



Gs+p{xi,x 2 ,x 3 ,x 4 ) = (0\Tip(x 1 )ip(x2)^(x 3 )if)(x 4 )\0) 

+ (0|T^(x 1 ) 7 5^2M^3)75^4)|0> , (4.1) 
g T (x 1 ,x 2 ,x 3 ,Xi) = {0\T$(x 1 )a^(x 2 )${x 3 )a' Jv iJ>(x 4 )\0) . (4.2) 

Using the general nonperturbative equation for the four-fermion function derived in the 
previous section, we can project out the equations for these S + P and T type four-point 
functions, which form closed set of linear integral equations. 

The four-fermion functions Qs+v and Qr are related to the functions Gs+p and Gt by 

Gs+p{Pl,P2,P3,Pi) = G S+ p(pi,P2,P3,P4) ~ ^G S +p(P3,P2,P1,P4) ~ 3Gt(P3, P2, Pi, Pi) (4.3) 
= G S +p{pi,P2,P3,P4) ~ -jGs+p{Pl,Pi,P3:P2) ~ SG T (pi, Pi, P3, P2) , (4-4) 

Gt(pi,P2,P3,P4) = G T (p 1 ,p 2 ,p 3 ,Pi) - ^G s+P (p 3 ,p 2 ,p 1 ,p4) + ^G T (p 3 ,p 2 ,p 1 ,p 4 ) (4.5) 
= G T (pi,p 2 ,p 3 ,p 4 ) - -G s+P (p 1 ,p 4 ,p 3 ,p 2 ) + -G T (pi,p 4 ,p 3 ,p 2 ) . (4.6) 

After the projection and reduction of the 7-algebra, we get for the s-channel S + P and T 
type four-fermion functions 

Gs+pijPi,P2,P3,Pi) = «e 2 ^V? / -^-\Z D ^(l) G s+p(Pi ~q,P2 + q,Ps,Pi) , (4.7) 

P1P2 J [ Z7T ) 

n ( \ ■ 2 ^P\Pl -Pi -JhcT f d 4 q 

G T {Pi,P2,P3,Pi) = ie ^2^2 J Jc^ d ^W g t\JPi ~ q,P2 + q,P3,P4) • (4.8) 

where D£(q) is the bare photon propagator. Remember that in Gs+p{pi,P2,P3,Pi) and 
Gt(Pi,P2,P3,P4), momentum conservation condition Pi+p 2 + p 3 + P4, = should be satisfied. 
Since the t-channel equations follow from these s-channel equations, it is enough to analyze 
these two s-channel equations, at least for the question whether or not there is a critical 
coupling. 

We also have to choose gauge. If we take the Feynman gauge, then the equation Eq. ( |4.8| ) 
gives immediately Gt(pi,P2,P3,P4) = 0, whereas the equation for Gs+p, Eq. ( |4.7|) , leads 
automatically to Gs+p(pi,P2,P3,Pi) = in the Yennie gauge. These results are obviously 
artifacts of the specific gauges, and one should choose the "correct" gauge in order to get 
reasonable results. It is known that the quenched ladder approximation for the fermion 
Schwinger-Dyson equation, which also follows from our effective action, is most reliable 
in the Landau gauge, a = 0. Therefore, we analyze the equations for these four-fermion 
functions in the Landau gauge, and use the photon propagator 

D,M) = ^ L, - q ~f) ■ (4-9) 
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A. G s+ P 



We start by looking at the equation for Gs+p, Eq. (J4.7|), which becomes in Euclidean 
space, after a Wick rotation, 



2 Pi - Pa /" d 4 g 1 



Gs+p(pi,P2,P3,P4) = ~3e —2-3- / TTT^I— G S+P(P1 - q,P2 + g,P3,P4) • 

P1P2 J \ A71 ) 1 

We redefine 

&s+p(Pi,P2,Ps,Pa) = G s+P (pi,p2 -P1,P3,P4) , 
which allows us to rewrite this equation into 



(4.10) 



(4.11) 



G s +p{pi,P2,P3,Pa) = -3e : 



•Pi ■ {P2 ~Pi) f d A q 1 



PHP2-P1) 2 J (2vr) 4 (q -pi) 2 



Gs+p{q,P2,P3,Pi) ■ (4.12) 



In order to solve this equation, we first expand G s+p(pi, P2, P3; Pi) m terms of p% 

00 

G'5+P(P1,P2,P3,P4) = X)(P2)w • • • (P2) Mn ^ 1 - At ' l (pi,p 2 r,P3,P4) , 



(4.13) 



n=0 



where we use the notation pf = pf/ypf for the angular dependence, and p ir = yjsf for 

the absolute value of the momenta. Note also that G^ 1 -" /in (pi,p 2r ,p3,p4) is symmetric in 
the indices pi, . . . , /x n . Inserting this in Eq. ( |4.12j ) shows immediately that the coefficients 
G^ 1 "^ n {pi,p2r,P3,Pi) have to satisfy 



PlG (Pl,P2r,P3,P4) = 3e 2 ' 



(2tt) 4 (g-p a ) 



-2pi 1 PiGo(pi,P2r,P3,P4) +p 4 Gf (pi,P2r,P3,Pi) = 3e 



2 G o(q,P2r,P3,P4) 
2 / A 



(27r)Mg-Pi) 2 

- pt'Goiq, p 2r , p 3 , p 4 ) + p 2 (5f (g, p 2r , Ps, P4) 



(4.14) 



(4.15) 



and 



^>?G^"( Pl ,p 2n p 3 ,P4) - 2pf p 2 G^-r(Pl,P2r,P3,P4) +p 4 ^ 1 -''^(pi,P2r,P3,P4) 

+ all symmetric permutations in . . . p n ) 
d A q 1 



3e 



pf <W" (g, p 2 r, P3, Pa) + p\G^ (g, p 2r , p 3 , p 4 ) 



(2vr)4(g-Pi) 2 
+ all symmetric permutations in {jL\ . . . /i n 

with n = 2, 3, 4, . . .. The homogeneous part of above equation satisfies 



p 2 iG^(p 1 ,p 2r ,p 3 ,p 4 ) =3e 2 



(2 7 r) 4 (g-p 1 ) 



r G^-^(g,p 2r ,p 3 ,P4). 



(4.16) 



(4.17) 
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In order to prove that the only solution of the equations Eqs. ( |4.14j )- (|4.16|) is the trivial 
solution, it is enough to prove that the only solution of the homogeneous equation Eq. Q4.171) 
is the trivial solution G^ 1 '"^ n {p 1 ,p2 r ,P3,Pi) = 0, for all n = 0,1,2,.... Note that only a 
solution which is symmetric in the Lorentz indices . . . fi n will contribute to the four-point 
function Gs+p ; from Eq. (|4.13|) it follows immediately that any nonsymmetric part does 
not contribute. 

The angular integration in Eq. ( J4.17[ ) can be performed analytically, using an expansion 
in a kind of generalized Chebyshev polynomials 



jj-fMl...fM2n 



"111 



($) = E 



-1*2 



kn2n—2k 



(2n - k)\ 



, (| k\(2n-2k)\(2n)\ 
+ all symmetric permutations in (/ix . . . /i 2n ) 



gMl qV2n-2k fit 1 2n-2k + ll J -2n-2k + 2 ( ^M2n-lA t 2n 



(4.18) 



E 



-l) k 2 2n+1 - 2k (2n + l-k)\ 



, n fc!(2n + l-2A;)!(2n + l)! 
+ all symmetric permutations in . . . fi2n+i 



q/J-l gM2n + l-2fc^M2n + 2-2fcM2n+3-2fc ^M2nM2n+l 

(4.19) 



where q 11 = / \fq 1 ~. These polynomials satisfy the identity 



da 



U^ n (q) 

' {q-Pi) 2 



n + l)p lr q r 



min(pi r , q r ) 
max(pi r , q r ) 



n+l 



(4.20) 



which is shown in Appendix |A], where we have given more details on these polynomials. 
Assuming that our four-point function is smooth enough (or more precisely, has a Taylor- 
expansion), we expandQthe unknown functions G^ 1 '" 1 ^™ (pi, p% r , p$, p±) in terms of these gen- 
eralized Chebyshev polynomials 



(4.21) 



1=0 



After this decomposition, which allows us to perform the angular integration, Eq. ( [4.17] ) has 
to hold for each function G™ separately, which leads to an equation for each coefficient 



pIgc::; 



(Plr,P2r,P3,P4.) 



3e 2 
16tt 2 



+ 



,p\ 



dq 2 Gf ul " M "(q r ,p 2r ,P3,P4)^ | • 



(4.22) 



For convenience we will drop all irrelevant indices and the spectator momenta in the 
discussion of this equation, defining g(p) as a scalar function of the absolute value of p\ 



2 This expansion can always be made in a unique way, and can be thought of as an expansion in 
p^ 1 . . .p^ m first, followed by an expansion of p^ 1 . . .pi™ in terms of U^ 1 " M '(pi) using Eqs. ( |A39 ) 
and ( A4C| ), see the appendix. 
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G fX.X (P lr > P 2r > ' = 9 i,P) 



(4.23) 



so the equation we have to consider becomes 

3e 2 



v g(p) 



16vr 2 



dq 2 g(q)+ r dq 2 ^- 2 g{ q )\ . 

Jp 2 q j 



(4.24) 



This integral equation is exactly the same as the equation for the dynamical mass function 
in quenched ladder QED, which has been studied extensively Ji4Hi~6|. It can be transformed 
into the following differential equation 



d d 
dp 2 dp* 



■p 2 g(p) 



3e 2 g(p) 
16tt 2 p 2 



(4.25) 



There are two independent solutions for the above differential equation, say (p 2 ) bi with bi 
satisfying 



(4.26) 



From this equation, it follows immediately that there is a critical coupling a c = tt/3, below 
which both 6, are real, and above which they are complex. 
So the general solution is 



9(p) = H °iP 



2b, 



(4.27) 



i=l,2 



in terms of the two independent solutions of the differential equation, with c, some arbitrary 
functions, depending on the spectator momenta pz r , P3, and p^ in the original equation 



Eq. ( |4.24| ), which needs to be determined. This general solution is constrained by the 
integral equation, and therefore, we substitute this solution back in the original integral 
equation, which gives us 



4 = 1,2 



c 26i+2 



(bi + I! 



+ p A 



K 2 bi 







(4.28) 



This boundary condition can only be satisfied above the critical coupling, if the exponents 
hi are complex; below the critical coupling, both bi are real numbers, and it is impossible 
to satisfy this boundary condition. Assuming there is a nontrivial solution, Eq. (|4.28|) is 
equivalent to the set of equations 



£1 

Cl 
C2 



6l + l 

~b 2 + l 
62 



.2(62-61) 



(4.29) 
(4.30) 



Below the critical coupling, this set of equations only has a solution if e > A, which means the 
infrared cutoff bigger than the ultraviolet cutoff. This is obviously not a physical solution, 
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and therefore the only solution of Eq. Q4.25|) is the trivial solution. That means that the 



only solution of the integral equation Eq ( 4.22Q is the trivial solution 



G ?7 1 ...u"(Plr,P2r,P3,P4) = . (4.31) 

Now going back to the original equation for the S + P type four-fermion function, using the 
decomposition Eqs. ( j4.21| ) and ( |4.13| ), and definition Eq. ( 4.11 ), we find that only the trivial 
solution 

Gs+p(PuP2,P3,Pa) = , (4.32) 

is allowed. 

Above the critical coupling, a > tt/3, both hi have a nonzero imaginary part, which 
implies oscillating solutions. Instead of using the complex exponents bi, we write the general 
solution of the differential equation Eq. (|4.25|) as 



g(p) = * sin (J *Z - 1 In 2) + Si cos (J^L - 1 In f) , (4.33) 



P \\ 71 AJ p \\j 71 A 

where we have used the ultraviolet cut-off A to set the mass scale. Substituting this solution 
back in original integral equation, we finally find a condition for A/e 

(2n + ^)tt + 5 n integer , (4.34) 




(4.35) 

Now one can always construct nontrivial solutions, depending on this ratio e/A. We assume 
that e, the infrared cutoff, is related to the energy scale of the dynamically generated con- 
densate through nonlinear effects, just as in the linearized quenched ladder approximation 
for the fermion Schwinger-Dyson equation, where the infrared cut-off sets the scale for the 
infrared mass function, e ~ m(0). Thus we allow the infrared cut-off to be a function of the 
ultraviolet cut-off, and in general, the lowest value of n for which e < A will correspond to 
the ground state solution; higher values of n correspond to oscillating solutions, i.e. excited 
states with higher energy. 

So there is a critical coupling a c = tt/3, above which one can find a nontrivial solution 
for the homogeneous equation, which is also the zeroth-order solution for G(pi,p 2 ,P3,P4) if 
one expand it in terms of p 2 - The only special point is that the ratio e/A must satisfy some 
constraint. In the limit that the coupling goes towards the critical coupling, a j a c , this 
ratio goes to infinity, and we can remove the ultraviolet and infrared cut-offs. 

Note that this critical coupling coincides with the critical coupling for dynamical mass 
generation for the fermions in the quenched ladder approximation. This is no surprise, since 
we noted already that we arrived at the same integral equation as that for the mass function, 
Eq ( f4.24|) . Looking back at the original equation for Gs+p, Eq. ( |4.10|) , we can understand 



this, because for p\ = p 2 this equation is exactly the same as the linearized quenched ladder 
approximation for the fermion propagator. 
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B. Gt 



Now we calculate Gt, Eq. ( |4.8[ ), which becomes in Euclidean space 
2 4 PiP2 -Pi -P2^ u f d A q q^q v 



G T (pi,P2,Pa,Pi) = e 



(2vr) 4 g 4 



Gt(pi ~q,P2 + g, Pa, P4) • (4-36) 



We can solve this equation in the same way as the S + P type equation, redefining 

Gt(pi,P2,P3,Pa) = G t (pi,P2-Pi,Ps,Pa) , 
which allows us to rewrite this equation into 



(4.37) 



G T (pi,P2,P3,P4) = e 



2 ^pt{p 2 - pi)" - pi ■ (p 2 - pi)g^ x 
3pi(P2 -Pi) 2 
d 4 g (g-pi)^(g-pi) v * , . 
(2vr) 4 (g-pi) 4 

Similarly, we expand Gt(px, p 2 , p 3 , Pa) in terms of p 2 

oo 

G t (Pi,P2,P3,Pa) = $3(Pa)i*i • • • ^^^'"""(PijParjPa,^) • 

n=0 



(4.38) 



(4.39) 



Inserting this in Eq. ( |4.38|) leads to a set of equations for the functions G f £ 1 -' tn (pi, p 2r , P3, Pa) , 
and in order to show that the only solution for Eq. ( 4.38| ) is the trivial solution, it is enough 
to prove that the only solution of the equation 



p 4 1 G^-^( Pl ,p 2r ,p 3 ,p i ) 
= e 2 (-4p?K + S^pl) 



d 4 q {% -Pi»){q v ~Pn 
(2vr) 4 {q-piY 



-G£ 1 - >ln (q,P2r,P3,P4) 



(4.40) 



is the trivial solution G% 1 "' fln (pi, p 2r , P3, Pa) = 0, for all n = 0, 1, 2, . . ., just as before. 
The generalized Chebyshev polynomials satisfy the identity 



(q-pi) 



SnO + U^iPl) 



71 + 1 



-q 2 - (pi - q 2 )p 



d 



lr 



lr 



PlrQV 



min(pi r , q r 
max(pi r , q r 



n+1 



(4.41) 



so if we expand the unknown functions G^ 1 "^ 71 (pi, p 2r , p 3 , p^) in terms of these polynomials 

oo 

G^( P l,P2r,P3,PA) =^Gf^(pir,P2r,P3,P±)Ur" Ul (Pl) , (4-42) 



1=0 



we can perform the angular integrals in Eq. ( J4.4U ) analytically, see Appendix This gives 
an equation for each of the coefficients separately 
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piG?:::: 



(Plr,P2r,P3,P4) 



487T 2 
r A 2 



d V 2G ?vi..M, (<lr,P2r,P3, P^tf 



5m + 



'1 + 2 l + 3q 2 



l + l l + lpi 



+ L dq 2 G?T'Z{(lr,P2r,P3,P±)q 2 



p'i 



$10 + 



I V \ I -I 
1 + 



l + lq 2 l + l 



PVr 
v Qr 



<h 

I 



(4.43) 



Again, for convenience we will drop all irrelevant indices and the spectator momenta in the 
discussion of this equation, defining g{p) as a scalar function of the absolute value of pi 



G fvZ^ (Plr , Plr , P3 , Pi) = 9 (Pi ) 



(4.44) 



so the equation we have to consider becomes 



p 4 g(p) 



48vr 2 

fA 2 



dq q g(q) 



no 



'1 + 2 l + 3q 2 



+ I dq 2 q 2 g{q) 



$10 + 



\l + l l + lp 2 
I p 2 I — 1 \ (i 



1 



+ 



l + lq 2 l + l 



(4.45) 



Since the I = and I ^ cases are different, we discuss them separately. 

For / = 0,the integral equation can be transformed into the following differential equation 



d d 
dp 2 dp r " 



p 6 g(p) 



e 2 

p 2 g(p) . 



167T 



(4.46) 



There are two independent solutions for the above differential equation, say (p 2 ) b 2 with b 
satisfying 



-1± Jl- 



a 



71 



(4.47) 



This means that the general solution of the differential equation can be written as 

9(p) = £ Q P 2iH ~ 2 , (4-48) 



i=l,2 



with Ci some arbitrary functions, depending on the spectator momenta p2 r , P3, and p^ in the 
original equation. We substitute this solution back in the original integral equation in order 
to get the correct boundary conditions, which gives us 



E « 

t=l,2 



h, - 1 



P 2 h 







(4.49) 



The critical coupling for the nontrivial solution is a° c = it, as can be seen from Eq. fl4.47|) , 
and just as for the S + P type equation, below the critical coupling this constraint allows 
only the trivial solution 



g(p) = 



(4.50) 
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Above this critical coupling, were the roots bi become complex, we can construct a nontrivial 
solution similar as before. 

For I ^ 0, the integral equation can be transformed into the following differential equation 



d d iP 2 Y +3 ^^(P 2 ) 



dp 2 dp 2 



dp 2 dp 2 



2 g(p) 



48vr 2 



d d 
dp 2 dp' 



(p 2 )^ 3 g(p) - 3(1 + 2)J- 2 (p 2 )^ 2 g(p) + 1(1 + 2)(p 2 )i +1 g(p) 



(4.51) 



Now there are four independent solutions for the above differential equation, say (p 2 ) b ' ~~ 
with b' satisfying the equation 



a 



//*(&') = b'(b' + 1)(6' -l)(b'-l-l) + — [3b'(b' - I - 1) + 1(1 + 2)] = . (4.52) 

127T 

We can calculate some special values /"(&')> to see where the roots are located, and find for 
small values of the coupling 



f?{l 4 
f?(l + 



2 

ft (I) 

no) 
//*(--) 



pa 
Ji 



a 

127T 
1 

_ 4 

a 

127T 

a 

127T 

1 

~4 

a 

12tt 



3^ 



/ 2 + 



1(1 + 2) > 

Z 2 + 2/ + ^- 
4 

1(1-1) > 
/(/ + 2) > 



a 1 12 71 9 



< 



(4.53) 



4 3tt 
(/ + 2)(/ + 3) > 



< 



From this it follows immediately that at small values of the coupling all four roots of the 
equation are real: two of them are positive and the other two are negative. In principle, 
there might be a critical coupling a c , where two roots become degenerate; such a critical 
coupling satisfies the equations 



/f(&') = o 
f[ a (b') = 0. 



(4.54) 
(4.55) 



This critical coupling will in general depend on the parameter I, and from Eq. (|4.53|) we can 
see immediately that a lower bound for these critical coupling is 



a\ > §71-, 
a 2 > i^vr 

"c — 53 ' 

a 3 > i^Tr 

c — 87 ' 



(4.56) 
(4.57) 
(4.58) 



A more detailed calculations shows that the critical coupling for I — 1 is indeed a\ = |tt. 
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Again, the general solution can be written as 

9(P)= E c lP 2b '^-\ (4.59) 



1=1,2,3,4 



with Q depending on the spectator momenta p2 r , P3, and P4 in the original equation 
Eq. (|4.43|) , and constrained by the integral equation. Therefore, we substitute this solu- 
tion back in the integral equation, which gives us 



4 



E C * 



0, (4.60) 



~l + 2 6^ 1 + 3 6 2 M+ 2 / pi+2 A 26<-2i-2 j _ j ^ A 2 ^~ 2; 

T+~\Ju l + / + lp ; + 2 (^ + l) ~ / + 1 ^ - / - 1 + / + 1 
which is equivalent to 

4 2b' 

£*-^ = o, ( 4 - 61 ) 

j=l °i 
4 e 2b^+2 

gc^^O, (4.62) 

4 A2M-2Z-2 

i=l U i 1 1 
4 7-1 a 26'-2Z 

We also have the additional constraint that the infrared cut-off is smaller than the ultraviolet 
cut-off. This set of equations has no nontrivial solution below the critical coupling, with all 
the roots b[ real, but above the critical coupling, some of b\ develop an imaginary part which 
means that the solution starts oscillating. Just as in the case of G$+p, these oscillations 
allow for nontrivial solutions of the integral equation, with e < A. 

Note, that the critical coupling for Gt depends on /, and all of them are bigger than that 
for Gs+p which is |; the smallest one, corresponding to I = is oP c = it. This means that 
one needs a stronger coupling to form a nontrivial T type function. If the coupling is in the 
region | < a < n, only Gs+p will have a nontrivial solution, while Gt still remains zero; 
at a = 7r this Gt might become nonzero as well. However, one should keep in mind that 
these Gs+pand Gt type function are not identical to the four-fermion functions Gs+p and 
Qt itself, see Eqs. ( |4.4p -(|OD. So as soon as either Gs+p (or Gt) becomes nonzero, both 
the S + P and the T type nonperturbative four-point functions Gs+p and Gt will become 
nonzero. 



V. DISCUSSION 

We have shown that, in leading order, both the S + P and the T type four-fermion 
functions and its condensates are zero below the critical coupling a c = | in the massless 



The reason to put the factor I — 1 in Eq. ( 4.64| ) is to keep the equations finite in the case of I = 1, 



in which case b' = 1 is a solution. We find in this case that in fact ^rrj- 
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fermion phase of a U(l) gauge theory. Above this critical coupling, there can be a nontriv- 
ial solution for these nonperturbative chiral-symmetry breaking four-point functions. This 
result is found using an expansion of the effective action, assuming that the other nonper- 
turbative four-fermion functions (which are not parity or Lorentz invariant) are zero. It 
agrees with the result found by Holdom and Triantaphyllou [fj], who used a truncation of 
the Schwinger-Dyson equation for the four-fermion functions. 

Our result shows that the critical coupling for chiral symmetry breaking by four-fermion 
condensates is the same as that for dynamical mass generation. Therefore, in a strongly 
coupled abelian gauge theory, beyond the critical coupling, the chiral symmetry will be 
broken through the two-fermion condensate since this is energetically favored. This means 
that in this theory, there is no hierarchy of chiral symmetry breaking, first by the four- 
fermion condensates, while the fermions remain massless, and another phase transition, at 
a lower energy, where the fermions acquire a mass. 

The generalization to the Nf flavor case is straightforward, and does not change our 
detailed calculations. The only thing we need to do is to label every Gi function with flavor 
indices. There are no flavor mixing effects in the equations, and the flavor symmetry is not 
believed to be broken dynamically, so the Green's functions for the different flavor indices 
decouple and the form of the equations is completely same as the ones we discussed. This 
means that in a Ui(Nf) <g> Un(Nf) gauge theory it seems to be impossible to use four- 
fermion functions or its condensates to substitute the two-fermion condensate as the chiral 
order parameter. 

In a more complicated theory like a chiral gauge theory, as suggested in ||, there exists 
in principle the possibility that chiral symmetry is broken through the forming of four- 
fermion condensates while the fermions remain massless. For such a scenario to be true, 
there should be some additional symmetry to prevent the fermions from becoming massive. 
Whether such a hierarchy of chiral symmetry breaking can indeed be realized in a chiral 
gauge theory, remains an open question; we only showed that in a U(l) gauge theory the 
critical coupling for the four-fermion condensates is the same as that for the two-fermion 
condensates. 

A similar result has been obtained using a truncation of the Schwinger-Dyson equation 
for the four-fermion functions 0. However, the problem with the set of Schwinger-Dyson 
equations is that it forms an infinite set of coupled equations, and it is not at all obvious 
how to truncate this set of equations in a consistent way. Using an approximation which is 
similar to the quenched ladder approximation for the fermion propagator Schwinger-Dyson 
equation, Holdom and Triantaphyllou have shown numerically that the critical coupling 
for chiral symmetry breaking by four-fermion condensates is almost the same as that for 
dynamical mass generation. It is remarkable that both their calculation and our analytical 
calculation give the same result (the same critical coupling), although the set of Feynman 
diagrams included in their and our truncation is not the same. This indicates that the 
essential diagrams are included in both schemes, and that the exact details of the truncation 
are less important. Of course, this should be confirmed by a study of the next-to-leading 
order effective potential, or a more accurate truncation scheme for the Schwinger-Dyson 
equation. We should also consider alternative formulations for the usual effective action, 
since there are different ways to introduce the auxiliary fields, which might lead to a different 
result. 
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APPENDIX A: ANGULAR INTEGRATION 

In this appendix we discuss the angular integration. Consider the case with three inde- 
pendent momenta, and choose the axes (in Euclidean space) to be 

Pi = Pi^ 

P2 = Pz(ni cos 6 P + n 2 sin 9 P ) (Al) 
q = q{n\ cos 9 qi + n 2 sin 9 qi cos 9 Q2 + n 3 sin 9 qi sin 9 q2 cos 4> q + sin 9 qi sin 9 q2 sin <p q ) , 

where we use the notation p for a four- dimensional vector and p for its absolute value. The 
angular integration we are interested in is 



q {q-pi) 2 



1 r a a r jo ■ 2n - a ^2 ( cos d P cos Ogi + sin d P sin d n cos 9 12) 



— I d9 n , I d9 n ~ sin 6L sin 9 n 



tt Jo qi Jo 92 91 92 q 2 + p\ - 2q Pl cos 9 qi 

q n p n . ^ (t 2 e - ie " + e w p) n+1 - (t 2 e ie f + e- i6 ") n+1 



<f dt(l~t 2 )t 

J unit circle 



2 n+3 7r(n + 1) sin 9 V Junit circle (q 2 + pi)t - qpx(\ + t 2 ) 

= 4,0 + 4,1 , (A2) 

where the last integration is performed along the unit circle in the complex t plane. There 
are two poles inside the circle, one is zero and the other is t±, which is defined as 



f± _ Pi + I 2 ± V (Pi + I 2 ) 2 ~ 4 Pig 2 _ f min[pi, g] max[pi, q] \ 
2piq \ max[pi, q) ' min[pi, q] J 

Their contributions to the integration are denoted by 4,0 and I n i respectively. We further 
find 



r n — l i 

rfh—ljyU ( L 2 J 

J ".o = - Ora+1 /, Jt -r- a E C k n+1 sm[{n+l-2k)9 p ] 



2 n+1 (n + l)pi sin# p 



' p\ + q 2 ^ 



fc=0,n>0 

rn + li 



£4-1- 



■2k 



' p\ + q 2 ^ 



fe=[f 



^ 2 + g 2N 
2pig 



E Cn+i sin[(n + 1 - 2k)9 p ]U n+1 „ 2k [ V -^- ) \ , (A4) 



„n—ln n+1 

/ n x = ^! p V- C k r n-l+2k in r (n + x _ 2jfe jfl 1 (A5) 

2™+ 1 (n + l)pisin0 p ^ n+1 ~ LV 1 pJ V ; 



where we have used the standard Chebyshev expansion JTy] 
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' p\ + Q 2 

1 -P^f t + t 2 ~ n \ 2 Pl q , 

PlQ 

Using another property of the Chebyshev functions 

sin[(n + 1)0] = U n (cos6) sin6 , 
the result of Eq. (|A2|) can be further expressed as 



f dn (q-P2) 2n+1 _ 9 2 >f +1 n ^ C n + l- k rj (cos0 )x Pl + q 2 \ 

J q (q-Pi) 2 ~ ^wrwi h 2n+2 c/2fe - i(cos ^ )x ^ 1 [-^r) 

where 

X n (x) = U n+1 (x) - U n ^(x) + (x- Vx^l) n+1 -(x + Vx^l) n+1 . 
To calculate X n (x), we consider the quantity 

oo oo 

J2 t n X n (x) = J2t n [Un+l(x) - £/„-i(x) + (X- Vx^l) n+1 -(x+ Vx^ 



n=0 n=0 
1 

~~ t 



1 - t 2 - 2t^x 2 - 1 

-1 + 



l-2xt + t 2 



x — vx 



Vx^T 



l-t(x- y/x 2 - 1) 



= 2£(x- Vx^l) n+1 t n , 

n=0 

where Eq. ( |A6| ) is used. Comparing both sides of the above equation, we find 



X n {x) = 2{x - Vx^l) n+1 . 
Now use the definition of the Chebyshev polynomials [IT 



fc =o k\{n-2k)\ 



to obtain 



{q-pir Piq t^o \ v\P2 J 

£W(*g) (-1)" " (2n+l) n ( Pi [ P2 \ 2k+2 

(q - p x ) 2 p x q ^ \ p x p 2 ) 
in which ajj. 2 "^ and a^ n+1 ^ are defined as 
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a 



(2n) 



= EH)" 



m=k 



a 



(2ra+l) 



E(-i) ? 



(n + m)! 
(n — m)!(m — k)\(m + 1 + fc)! ' 
(n + m + 1)! 



(n — m)!(m — k)l(m + 2 + fc)! 
Finally, to calculate the coefficients a, we introduce two functions 

(to + m + 2A;)! 



F n>k {x)-j2{ i r (n _ m)!m!(m+1+2A , )! 

n,feW ~^ J (TO-m)!m!(m + 2 + 2A;)! 



so we get 



(2n) 
(2n+l) 



For to < 1 we find 



Fo,k( x ) 



l) fc F n „ fc)fc (l) 
v — l) fc F n _fc i fe(l) 

1 



2fc + 1 ' 

^ (x) = 2TT^' 

F Xifc (x) = 1 - x, 
Fi t k(x) = l-x, 



whereas for to > 1 



Fn.k\X] .X 



-2k- 



_ x d 



n-l 



—y k+n {i-xY\ 

I d n ~ x 



r?! 



F n , fc (x) = -x- 2fc - 2 -— T [x 2fc+ " +1 (l -x)»] 



so the total result can be written as 

F n>fe (l) 
F n , k (l) 



(2n) 



2fc + 1 ' 

<5 n 
2TT2 ' 



a 



(2n+l) 



2to + 1 

(-1)" 
2to + 2 



■6. 



nk 



$nk 



Substituting this result into Eqs. ( |A14j ) and ( |A15 ), we finally get the result 



'(q-pi) 2 (n + l)piq 



min[px, q] 
max[pi, q\ 



n+l 



U n (Pi - Pa) 
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where p = p/p denotes the angular dependence. 

Next, we introduce a "generalized Chebyshev function" with Lorentz indices fii . . . fi n by 
extracting the p dependence in U n (q ■ p) 



U n {q-p)=p^...pjj^{q), 
which defines the generalized Chebyshev functions U^ 1 '" tln as 



(A33) 



'2ll 



_ ™ (-l) k 2 2n ~ 2k (2n - k)\ 
~ ^ k\(2n - 2k)\(2n)\ 

+ all symmetric permutations in (fii . . . /i 2n 

(-\) k 2 2n+1 - 2k (2n+l-k)\ 



, „ fc!(2n + l-2A;)!(2n + l)! 
+ all symmetric permutations in (/ii . . . /i2n+i 



gVl gf J -2n-2k fi[ J -2n-2k + ll J -2n-2k + 2 ( ^M2n-lA t 2n 

(A34) 

gMl gA t 2n + l-2fc ( ^M2n + 2-2fcA t 2n+3-2fc ^M2nM2n + l 

(A35) 



Since in Eq. ( |A32| ), P2 is an arbitrary momentum variable, we can use the expansion 
Eq. ( |A33| ) to find 



da 



' (q-pi) 2 



(n + l)piq 



min[pi, q] 
max[pi, q] 



n+l 



(A36) 



This formula allows us to perform all of the angular integration which are needed, rest 
us to show that we can indeed expand our four-point functions in terms of these generalized 
Chebyshev functions. The Chebyshev polynomials U n (x) form a complete set of indepen- 
dent functions, so from Eq. ( |A33[ ), it follows immediately that our generalized Chebyshev 
functions are also a complete set for functions which are totally symmetric in their Lorentz 
indices. The explicit expression of this statement can be seen from following identities, 
noting that 



(P • Q) 2n = E Cn,kU 2k (p ■ q) 



k=0 



Pn ■ • -2W E c ntk U^^(q)5^^ . . . ) 



(A37) 



k=0 



(p-q) 



2n+l 



E C n,k U 2k+l(p ■ q) 



k=0 



n n r ' fj^-^k+l ( fi\ S^2k + 2^2k + 3 ,W2nM2n4 

p^i ■ ■ ■ Pfj.2n+i °n,k u 2k+i yq)° ■ ■ ■ ° 

k=0 



(A38) 



which follows directly from Eq. ( |A33| ) and the fact the Chebyshev polynomials form a com- 
plete set. Now due to the fact that p is arbitrary, the Eqs. ( |A37| ) and ( |A38| ) are equivalent 
to 



\ -y Cn,k 



Jjill—illk ^fiV2k+lV2k+2 ^A12n-lM2n 
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+ all symmetric permutations in (fj,i . . . ji2n) 



-7/^1 



7M2n + l 



E 



u 2k+l 



^£M2fc+2M2fc+3 £M2n/i2n+l 



,=o(2n+l)! 
+ all symmetric permutations in (//! . . . //2n+i) 



(A39) 



(A40) 



So all g m . . . g Mn can be expanded in terms of our generalized Chebyshev polynomials. Since 
any smooth function of q can always be expanded in terms of g w . . . q^ n times a function 
which depends only on the absolute value q, this means that we can always expand our 
four-point functions in terms of these Chebyshev functions. 
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FIGURES 

FIG. 1. Our truncated equation for the nonperturbative four-fermion function (diagrammati- 
cally): (a) G(pi,P2,P3,P4) decomposed into the s- and ^-channel components G(pi,p2,P3,P4); (b) 
the equation for the s-channel component; and (c) the equation for the t-channel. 
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Figure 1(c) 



